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Abstract 

A multitype Dawson- Watanabe process is conditioned, in subcriti- 
cal and critical cases, on non-extinction in the remote future. On every 
finite time interval, its distribution is absolutely continuous with re- 
spect to the law of the unconditioned process. A martingale problem 
characterization is also given. Several results on the long time behav- 
ior of the conditioned mass process — the conditioned multitype Feller 
branching diffusion — are then proved. The general case is considered 
first, where the mutation matrix which models the interaction between 
the types, is irreducible. Several two-type models with decomposable 
mutation matrices are also analyzed. 
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Introduction 



The paper focuses on some conditioning of the measure- valued process called 
multitype Dawson- Watanabe (MDW) process, and on its mass process, the 
well-known multitype Feller (MF) diffusion. We consider the critical and 
subcritical cases, in which, for any finite initial condition, the MF diffusion 
vanishes in finite time, that is the MDW process dies out a.s. In these cases, 
it is interesting to condition the processes to stay alive forever - an event 
which we call remote survival, see the exact definition in ([3]). 
Such a study was initiated for the monotype Dawson- Watanabe process by 
A. Rouault and the second author in [2H] (see also [TU], [U] and [H] for the 
study of various aspects of conditioned monotype superprocesses). Their 
results were a generalization at the level of measure-valued processes of the 
pioneer work of Lamperti and Ney ([20], Section 2), who studied the same 
questions applied to Galton- Watson processes. 

We are interested here in the multitype setting which is much different 
from the monotype one. The mutation matrix D introduced in ([2]), which 
measures the quantitative interaction between types, will play a crucial role. 
We now briefly describe the contents of the paper. The model is precisely 
defined in the first section. In the second section we define the conditioned 
MDW process, express its law as a locally absolutely continuous measure 
with respect to the law of the unconditioned process, write explicitly the 
martingale problem it satisfies and give the form of its Laplace functional; 
all this in the case of an irreducible mutation matrix. Since D is irreducible, 
all the types communicate and conditioning by remote survival is equivalent 
to conditioning by the non-extinction of only one type (see Remark 12.51) . 
The third section is devoted to the long time behavior of the mass of the 
conditioned MDW process, which is then a conditioned MF diffusion. First 
the monotype case is analyzed (it was not considered in [28]), and then the 
irreducible multitype case. We also prove that both limits interchange: the 
long time limit and the conditioning by long time survival (see Theorem 13.71) . 
In the last section we treat the same questions as in Section [H] for various 
reducible 2-types models. Since D is decomposable, the two types can have 
very different behaviors, that also depend on the precise conditioning that is 
considered (see Section H?T1) . 

1 The model 

In this paper, we will assume for simplicity that the (physical) space is R. 
k is the number of types. Any /c-dimensional vector u G is denoted by 
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■ ■ ■ ; Uk). 1 will denote the vector (1; . . . ; 1) G R''. \\u\\ is the euclidean 
norm of u E R'' and {u,v) the scalar product between u and v in [R^. If 
M G R^, |m| is the vector in IR^ with coordinates \ui\, 1 < i < k. 
We will use the notations u > v (resp. u > v) when u and f are vectors or 
matrices such that u — v has positive (resp. non-negative) entries. 
Let Cb(R, R^) denote the space of [R'^-valued continuous bounded functions 
on R. By C;,([R, [R*^)+ we denote the set of non- negative elements oi CbiR^R^). 
M{R) is the set of finite positive measures on R, and M{RY the set of k- 
dimensional vectors of finite positive measures. 

The duality between measures and functions will be denoted by (■, •) : 
(z/, /) := J fdv if z/ G M(R) and / is defined on R, and in the vectorial 
case 

((i/i;...;z/fc),(/i;...;/fc)) fid^i = {{{i^i, fi); . . . ; {uk, fk)),l) 

1=1 ^ 

for u = (z/i; . . . ; z/fc) G M(R)'= and f = {f,;...; /,) G ^(R, R')- For any 
A G R'^, the constant function of Cb{R, R^) equal to A will be also denoted by 
A. 

A multitype Dawson- Watanabe process with mutation matrix D = 
i<i,j<k is a continuous M(R)'^-valued Markov process whose law P on 
the canonical space {Q := C(R+, M(R)'^), {Xt)t>o, {^t)t>o) has as transition 
Laplace functional 

V/ G a(R, R^)+, E(exp -(X„ /) | Xq = m) = exp -(m, UJ) (1) 

where Utf G Cfe(R, R'^)+, the so-called cumulant semigroup, is the unique 
solution of the non-linear PDE 

= AUJ + DUJ - ^{Ujr (2) 
Uof =/. 

Here, u Q v denotes the componentwise product (ujfi)i<i<fc of two k- 
dimensional vectors u and v and u®^ = u Q u. To avoid heavy notation, 
when no confusion is possible, we do not write differently column and row 
vectors when multiplied by a matrix. In particular, in the previous equation, 
Du actually stands for Du'. 

The MDW process arises as the diffusion limit of a sequence of particle sys- 
tems {^N^)x, where is an appropriate rescaled multitype branching 
Brownian particle system (see e.g. [13] and (TH], or [32] for the monotype 
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model): after an exponential lifetime with parameter each Brownian par- 
ticle splits or dies, in such a way that the number of offsprings of type j 
produced by a particle of type i has as (nonnegative) mean 5ij + j^dij and 
as second factorial moment c {5ij denotes the Kronecker function, equal to 1 
ii i = j and to otherwise). Therefore, the average number of offsprings of 
each particle is asymptotically one and the matrix D measures the (rescaled) 
discrepancy between the mean matrix and the identity matrix /, which cor- 
responds to the pure critical case of independent types. 
For general literature on DW processes we refer the reader e.g. to the lectures 
of D. Dawson [3J and E. Perkins [25] and the monographs [5] and [7]. 
Let us remark that we introduced a variance parameter c which is type- 
independent. In fact we could replace it by a vector c = (ci; ■ ■ ■ ; c^), where 
Cj corresponds to type i. If infi<j<fcCj > 0, then all the results of this paper 
are still true. We decided to take c independent of the type to simplify the 
notation. 

When the mutation matrix D = {dij)i<i j<k is not diagonal, it represents 
the interaction between the types, which justifies its name. Its non diagonal 
elements are non-negative. These matrices are sometimes called Metzler- 
Leontief matrices in financial mathematics (see [29] § 2.3 and the bibliography 
therein). Since there exists a positive constant a such that D + aJ > 0, it 
follows from Perron- Frobenius theory that D has a real eigenvalue /i such 
that no other eigenvalue of D has its real part exceeding /x. Moreover, the 
matrix D has a non negative right eigenvector associated to the eigenvalue 
II (see e.g. [H], Satz 3 § 13.3 or [29] Exercise 2.11). The cases /i < 0, /i = 
and /i > correspond respectively to a subcritical, critical and supercritical 
processes. 

In the present paper, we only consider the case /i < 0, in which the MDW 
dies out a.s. (see Jirina [T7]). 

2 (Sub) critical irreducible MDW process 
conditioned by remote survival. 

Let us recall the definition of irreducibility of a matrix. 

Definition 2.1 A square matrix D is called irreducible if there is no per- 
mutation matrix Q such that Q^^DQ is block triangular. 

In all this section and in the next one, the mutation matrix D is assumed 
to be irreducible. By Perron- Frobenius' theorem (see e.g. [2H] Theorem 1.5 
or [H], Satz 2 §13.2, based on [2j and [13]), the eigenspace associated to 
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the maximal real eigenvalue of D is one-dimensional. We will always de- 
note its generating right (resp. left) eigenvector by ^ (resp. by rj) with the 
normalization conventions (^,1) = 1 and = 1- All the coordinates of 

both vectors ^ and rj are positive. 



2.1 The conditioned process as a /i-process 

The natural way to define the law P* of the MDW process conditioned to 
never die out is by 

VE G ^i, r{B) := lim P{B \ {Xt+e, 1) > 0) (3) 
if this limit exists. 

The following Theorem 12.21 proves that P* is well-defined by ([3]) and is 
a probability measure on J^t absolutely continuous with respect to P. Fur- 
thermore, the density is a martingale, so that P* can be extended to \/t>o^ti 
defining a Doob /i-transform of P (see the seminal work [22] on /i-transforms 
and [21] for applications to monotype DW processes). 



Theorem 2.2 Let P he the distribution of a critical or suhcritical MDW 
process characterized by ([T]), with an irreducible mutation matrix D and ini- 
tial measure m G M(P)'^ \ {0}. Then, the limit in exists and defines a 
probability measure P* on \/t>o^t such that, for any t > 0, 



p*| =i^^e-Mt p| f4) 



where ^ E P!' is the unitary right eigenvector associated to the maximal real 
eigenvalue of D. 



Proof of Theorem 12.21 By definition, for B E J^t, 

E(lB(l-P((Xi+,,l) = 0|.Ft))) 



E(1b I {Xt+e, 1) > 0) 



P((Xi+e,l) = 0) 



For any time s > 0, := ((X^,!, 1); • • • ; 1)), the total mass at time 

s of the MDW process — a multitype Feller diffusion — is a continuous Re- 
valued process with initial value x = ((mi, 1); . . . ; (m^, 1)) characterized by 
its transition Laplace transform 

VA e Pi, E(e-("*'^) \xo = x) = e-("'"'\ (5) 
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Here, = (u^i; . . . ; u^fc) := UtX satisfies the non-linear differential system 



dt t ^\ tl (6) 



or componentwise 



du^ c 



Then, 



P((X„ 1) = 0) = lim E(e-<^-^> | Xq = m) = e-^^'""^^-- 

A^oo 



where A oo means that all coordinates of A go to +oo. Using the Markov 
property of the MDW process, one obtains 

E(i, I 1) > 0) = \_^,„„„,_„.^., m 



In the monotype case {k = 1), can be computed explicitly (see Section 
I3.1I) . but this is not possible in the multitype case. Nevertheless, one can 
obtain upper and lower bounds for u^. This is the goal of the following 
two lemmas, the proofs of which are postponed after the end of the proof of 
Theorem I2.2[ 

Lemma 2.3 Let = {u^^; . . .; u^^.) he the solution of 

(i) For any \eR\\ {0} and any t > 0, > 0. 

■ 

(ii) Let := sup — ^ and ^ := infj^j. Fort > and A G R^, 

- in the critical case (fi = 0) 

< ^ and therefore sup < — (8) 

- in the subcritical case (fi < 0) 

< ^ and therefore sup < (9) 
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fiii) Let B 



A 



U 



inf — ^ and ^ :=sup, fj. Then 
i<i<k ^ 



Vt > 0, AG R'l, 



1 + fB^ t 

1 + ^^^0 (1 - en 



if 12 = 
if jj, <0. 



(10) 



(iv) For any A G IR^ and t > 0, 



1 + |Co^ t 



,Dt ■ 



if 11 = 



1 + (1 



e^*A z//i<0. 



(11) 



The main difficulty in the multitype setting comes from the non- 
commutativity of matrices. For example (EI) can be expressed as ^ = 
[D+At)u^ where the matrix At is diagonal with z-th diagonal element 
However, since D and At do not commute, it is not possible to express in 
terms of the exponential of Jq{D + Ag) ds. The following lemma gives the 
main tool we use to solve this difficulty. 

Lemma 2.4 Assume that 1 1— > f{t) G IR a continuous function on [R+ and 
t Ut €i R'' is a difjerentiable function on IR_|_. Then 

'^>{D+f{t)I)ut, Vt>0 = 



ut > exp (^^ {D+f{s)I) ds^uo, Vt > 



For any I < i < k, applying ((51) with x 



e* where e* = 6ij, 
one easily deduces the existence of a limit in [0, oo] of when A 
Moreover, by Lemma [2.31 (ii) and (iii), for any t > 0, 



l<J<k, 
oo. 



< < lim < < +00 

A^oo 



where f{e) = 1/0 if = or /(0) = |yu|e^7(l - e^^) if /i < 0. Therefore 
lime^oo liniAc-^oo Uq=0 and, for sufficiently large 6*, 



1 



g-(xt,limA^oo «e) 



1 



< K 



for some constant that may depend on t but is independent of 6. Since 
E(Xt, 1) < cx) for any t > (see [I5] or [16]), Lebesgue's dominated con- 
vergence theorem can be applied to make a ffist-order expansion in 6' in ([71). 
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This yields that the density with respect to P of P conditioned on the non- 
extinction at time t + 6 on J-'t, converges in L^{P) when 6' oo to 

if this hmit exists. 

We will actually prove that 

1 e~^* 
lim sup II- — M^-— — ^11=0. (13) 



This will imply that the limits in 9 and in A can be exchanged in (]T2i) and 
thus 

lim E(1b I {Xt+e. 1) > 0) = e"^* cfl^^) = e"^* pfl^^^), 

completing the proof of Theorem 12. 2[ 
Subcritical case: < 

As a preliminary result, observe that, since D has nonnegative nondiagonal 

entries, there exists a > such that D + a/ > 0, and then exp(Dt) > 0. 

dii^ 

Since — < Du'^, we first remark by Lemma [2.41 (applied to —u^), that 



dt 



Vt > 0, < e^*A. 



Second, it follows from Lemma 12.31 (iv) that 



e^*A - < (^1 - (1 + ^C'o"(l - e'^*))^«"/^^e^*A 

< ^Co^(i_eA^*)e^*A. 

Therefore, since Cq = supj Aj/^,, there exists a constant i^' independent of A 
such that 

VA > 0, Vt > 0, e^*A - < ir||A||e^*A. (14) 

In particular, 

||A|| < ^7 -e^*A Vt > 0. 

" " - 2ir * - 2 

Third, it follows from Lemma 12.31 (ii) that there exists to such that 

Vt>to, VA>0, ||m,^|| <^. 
8 



Fourth, as a consequence of Perron-Frobenius' theorem, the exponential ma- 
trix e^* decreases hke e^* for t large in the following sense: as t — > oo, 

3^ > 0, e^* = e^*P + ©(e^'^-^)*) (15) 

where P := {^iVj)i<i,j<k (see ^29l| Theorem 2.7). Therefore, there exists 6*0 
such that 

Vt > 00, ^e'^^P < e^* < 2e''*P. 
Last, there exists a positive constant K' such that 

Vm,wGR^, (t;,PM) = (u,r7)(t;,0 > ^^(m, 1) (^^, 1) > 

Combining all the above inequalities, we get for any a G IR^, 6, A G IR^ \ {0} 
and for any 9 > 9q, 

(a,<+,) (a,e™0 



(&,e^(^+*)<) 

(a,|<^,-e^Xl) , Ke™<)(6,|<^.,+,-e^(^+*)<|) 



< 



(6,e^(^+*)^J (6,e^(^+*)02 
< Allan IK lie l(^^p,A)+ (^^p^^). 



<^e-*yil<ll (16) 

where the constants K may vary from line to line, but are independent of A 
and to- 

Now, let to(^) be an increasing function of 6 larger than such that 
^o(^) oo when 6 oo. By Lemma [2.31 (ii), —>■ when 9 oo, 

uniformly in A > 0. Then, by ffTUj) . uniformly in A > 0, 

lim ^^''^^oid)+e) _ ^^^t'oCe)) 



^--(&,eM^+*)P<(,)) 



e 
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which completes the proof of Theorem 12.21 in the case fi < 0. 
Critical case: /i = 

The above computation has to be slightly modified. Inequality ( fT^ becomes 

\uf - e^*A| < (^1 - (1 + |Co' t)~^"/^^e^*A 

<ir||A||te^*A. (17) 
Therefore, the right-hand side of (fT6l) has to be replaced by 

+ (18) 

Now, using Lemma 12.31 (iii) again, it suffices to choose a function to(^) in 

such a way that \img^oo(^ s,upx>o\\uf^(^e)\\ ~ ^- '^^^ ^'^^ complete the 

proof of Theorem 12.21 as above. □ 



Proof of Lemma l2.3l 

(i) First, observe that, by ([5]), > for any t > 0. Next, since D is 
nonnegative outside the diagonal, 

k 

Therefore, for any i such that Aj > 0, u^^ > for any t > 0. 
Let I := {i : Xi > 0} and J := {j : Xj = 0}. By the irreducibility of the 
matrix D, there exist i ^ I and j & J such that dji > 0. Therefore, for 
sufficiently small t > 0, 

1=1 

and thus u^j > for t > in a neighborhood of 0. Moreover, as long as 

> 0, for the same reason, u^ j cannot reach 0. 
Defining /' = I U {j} and J' = J \ {j}, there exists i' G /' and j' G J' 
such that dj'i' > 0. For sufficiently small e > 0, m^^, > and the previous 
argument shows that u^+tj' > for t > as long as u'^+^y > 0. Letting e go 
to yields that u^^, > for sufficiently small t > 0. 

Applying the same argument inductively shows that > for t > in a 
neighborhood of 0. Using (fTOll again, this property can be extended to all 
t > 0. 



(Hi 
dt 
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(ii) and (iii) As the supremum of finitely many continuously differen- 
tiable functions, t t— > is differentiable except at at most countably many 
points. Indeed, it is not differentiable at time t if and only if there exist two 
types i and j such that u^^i/^i = u^ j/^j and d{u^ J^i)/dt ^ d{u^j/^j) /dt. 
For fixed i and j, such points are necessarily isolated, and hence are at most 
denumerable. 

Fix a time t at which is differentiable and fix i such that u^^ = C^^i. 
Then 



dC^ ^ duf. 



dt dt l^'^'^^'^'J 2^ 



where the inequality comes from the fact that D is nonnegative outside of 
the diagonal and where the third line comes from the specific choice of the 
subscript i. Therefore, 

^<^^C,'-U^C,r- (20) 

Assume fi = 0. 

By Point (i), if A 7^ 0, > for any t > (the case A = is trivial). Then, 
for any t > 0, except at at most countably many points, 

dCt^/dt ^ c ^ 



Integrating this inequality between and t, we get 



C,'-C^ 2^ '-l + 23Lj(' 

The proof of the case /i < can be done by the same argument applied to 
t 1-^ e~^^C^. Inequalities (iii) are obtained in a similar way too. 
(iv) By definition of C^, ([6]) implies that 



dt 

Then, (iv) follows from (ii) and Lemma [2. 4[ □ 
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Proof of Lemma 12.41 Fix e > and let 

u'f^ := exp { [ {D + f{s)I)ds){uo - e). 
Jo 



Then 



dt dt 



Let to •= inf{t > : 3z G {1, . . . , k},Ut^i < uf}}. For any t < to, since D is 
nonnegative outside of the diagonal, 

e {1, . . . , k}, j^iut, - u^:}) > {du + fmut, - uf}). 

Since uq > ul^\ this implies that Ut — uf^ > for any t < to and thus 
to = +00. Letting e go to completes the proof of Lemma [2^ □ 

Remark 2.5 Since the limit in ( fT3l) is uniform in A, one can choose in par- 
ticular X = X\l < i < k, where A* = for j ^ i. Thus, for each type 

h 

hni — = 4 

which implies as in that, for B & J-'t, 

lim P{B I {Xt+e,i, 1) > 0) = lim P{B \ {Xt+e, 1) > 0) = P*(5). 

Therefore, Theorem 12.21 remains valid if the conditioning by the non- 
extinction of the whole population is replaced by the non-extinction of type 
i only. This property relies strongly on the irreducibility of the mutation 
matrix D. In Section |U we will show that it does not always hold true when 
D is reducible (see for example Theorem 14.11 or Theorem 14. 4p . ^ 

2.2 Laplace functional of P* and Martingale Problem 

To better understand the properties of P*, its Laplace functional provides a 
very useful tool. 

Theorem 2.6 P* is characterized by: V/ G Cfe(R, R^)+ 

r (exp -{X„ f) I Xo = m)= ^J^^pYl^^-^^^-^^^^^f) (21) 
inhere the semigroup Vtf is the unique solution of the PDF 

^ = AVJ + DVJ - cUJ © VJ, Vof = e (22) 
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Proof From Theorem 12.21 and ([T]) we get 

\{m,0 J 



-fit a 

-(m,Utf)_:±_ 



e=0 



e -f"^ /™ rr A d 



(m, 4) oe '^-'J 

Let Vtf := §-Ut{f + e^l^^- Then y,/ is solution of 



at 



£=0 



andT/o/ = f(/ + ^OL=o = ^- ° 

Comparing with the Laplace functional of P, the multiplicative term 
^"^m^i) appears in the Laplace functional of P*. In particular, the multi- 
type Feller diffusion Xt is characterized under P* by 

E*(exp -(xi, \)\x^ = x) = i^^e-^*e-(^'"'), A G (23) 

where := Vt\ satisfies the differential system 

-^ = Dv^~cufQv^, v',=t (24) 

The following theorem gives the martingale problem satisfied by the con- 
ditioned MDW process. This formulation also allows one to interpret P* 
as an unconditioned MDW process with immigration (see Remark 12.81 be- 
low). The term with Laplace functional e~^* that we just mentioned 
is another way to interpret this immigration. 

Theorem 2.7 P* is the unique solution of the following martingale problem: 
for all f G Ci{R, P'=)+, 

exp(-(Xt, /)) - exp(-(m, /)) 

+ £ {{Xs, (A + D)f) + c^-^jj/^ - ^{Xs, f^')) exp(-(X„ /)) ds 

(25) 

is a P* -local martingale. 
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Proof According to [15J (see also [B] for the monotype case), P is the unique 
solution of the following martingale problem : for any function F : M(IR)'^ 
R of the form /)) with G C2(R, R) and / G Cl{V, R'=)+, 

F{Xt) - F{Xo) - [ AF{Xs) ds is a P-local martingale. (26) 

Here the infinitesimal generator A is given by 

dF c 

AFim) = (m, (A + D)—) + -(m, d^F/dm^) 

dm 2 

.dF dF ^ c A, S^F, 

i=l * j=l J i=l * 

where we use the notation OF /dm = [OF / dmi)i<i<k and d'^Fjdm? = 
{d'^F/dm})i<i<k with 

dF , 1 

'x) := lim - (F(mi, . . . , mj + eS^, • • • , 'n^/c) — -F(7Ti)) , x G R. 



dmi e^o e 

Applying this to the time-dependent function 

F(s,m) := (m,e)e-'^^e-<"^'^> with / G ^^^(R,^'^ 

for which 

gF(g,m) ^ _/(a;)F(s,m) +ee-^^-<'"'^> 
am 



and 



d F{s, m 



dm^ 
one gets 

dF 

— (s,m)+^(F(s,-))(m) 



l(a;) = f®\x)F{s, m) - 2f{x) ^e''^^-^^^'^^, 



-(m, (A + Z})/)F + |(m, - c ^^;^®^^ F. 



Therefore, 



+ e)e--^((X„ (A + D)f) + c^^^^ - |(X„ /«^))e-<^-^)rf. 

is a P-local martingale, which implies that (!25l) is a P*-local martingale. 
The uniqueness of the solution P* to the martingale problem ( l25l) comes from 
the uniqueness of the solution of the martingale problem (126|) . □ 
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Remark 2.8 Due to the form of the martingale problem the probabil- 
ity measure P* can be interpreted as the law of a MDW process with inter- 
active immigration whose rate at time s, if conditioned by X^, is a random 
measure with Laplace functional exp — c-^^^^^^. Monotype DW processes 
with deterministic immigration rate were introduced by Dawson in [2]. The 
first interpretation of conditioned branching processes as branching processes 
with immigration goes back to Kawazu and Watanabe in [131, Example 2.1. 
See also [2H] and [lU] for further properties in the monotype case. <C> 



3 Long time behavior of conditioned multi- 
type Feller diffusions 

We are now interested in the long time behavior of the MDW process con- 
ditioned on non-extinction in the remote future. Unfortunately, because of 
the Laplacian term in (|22l) . there is no hope to obtain a limit of Xt under P* 
at the level of measure (however, see flO] for the long time behavior of crit- 
ical monotype conditioned Dawson- Watanabe processes with ergodic spatial 
motion). Therefore, we will restrict our attention to the R'^- valued multitype 
Feller diffusion Xt- As a first step in our study, we begin this section with 
the monotype case. 



3.1 Monotype case 

In this subsection, we first give asymptotic behavior of Xt under P* (Propo- 
sition [3]T]). This result is already known, but we give a proof that will be 
useful in the following section. We also give a new result about the exchange 
of limits (Proposition 13. 3p . 

Let us first introduce some notation for the monotype case. 
The matrix D is reduced to its eigenvalue fi, the vector ^ is equal to the 
number 1. Since we only consider the critical and subcritical cases, one has 
fi < 0. The law P* of the MDW process conditioned on non-extinction in the 
remote future is locally absolutely continuous with respect to P (monotype 
version of Theorem l2.2l already proved in [2H], Proposition 1). More precisely 



P*| {^t, 1) p| (27) 
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Furthermore the Laplace functional of P* satisfies (see PS], Theorem 3): 

r(exp-(X„/) I Xo = m) = i^^e-'^^e-^-^*/) = i^^e-^^'^*/) 

(m, 1) (m, 1) 

(28) 

where 

^ = AVtf-cUtfVtf, Vof = l. (29) 

The total mass process Xt = {Xf, 1) is a (sub) critical Feller branching 
diffusion under P. By ( l28l) its Laplace transform under P* is 

E*(exp -\xt \xo=x)= e-^'^'s A G P+, (30) 

with 

Recall that the cumulant satisfies 

^ = /^«*^-|K^)^ 4 = x. (31) 

This yields in the subcritical case the explicit formulas 



and = exp ( —c [ u^ds] = 

\ Jo J fl + ^A(l- 



A > 0, (32) 

(33) 



In the critical case (/i = 0) one obtains (see [20| Equation (2.14)) 



uf = ^ , and = = ^. (34) 



We are now ready to state the following asymptotic result. 
Proposition 3.1 

(a) In the critical case (fi = 0), the process Xt explodes in ¥* -probability when 
t oo, i.e. for any M > 0, 

lim P*{xt < M) = 0. 

t— >+oo 
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(b) In the subcritical case (fi < 0), 

lim P*(xt G ■) = r(2,^) 

where this notation means that Xt converges in P* - distribution to a 
Gamma distribution with parameters 2 and 2\fi\/c. 

One can find in [19] Theorem 4.2 a proof of this theorem for a more general 
model, based on a pathwise approach. We propose here a different proof, 
based on the behavior of the cumulant semigroup and moment properties, 
which will be useful in the sequel. 



Proof For /i = 0, by (!34|) . -u^ — > and — when t — oo for any 
A 7^ 0. This implies by (1301) the asymptotic explosion of Xt in P*-probability. 
Actually, the rate of explosion is also known: in [TU] Lemma 2.1, the 
authors have proved that y converges in distribution as t ^ oo to a 
Gamma-distribution. This can also be deduced from since u^^^ and 

v^^^ converge to and 1/(1 + cA/2)^ respectively, as t ^ oo. 

For /i < 0, by (l30l) . ( l32l) and ( l33l) . the process Xt has the same law as 
the sum of two independent random variables, the first one with distribution 
r(2, ^(j^i^pi) ) and the second one vanishing for t ^ oo. The conclusion is now 
clear. □ 



Remark 3.2 The presence of a Gamma-distribution in the above Proposi- 
tion is not surprising. 

• As we already mentioned it appears in the critical the limit law 
of xt/t [lO]. 

• It also goes along with the fact that these distributions are the equilib- 
rium distributions for subcritical Feller branching diffusions with con- 
stant immigration. (See pj, and Lemma 6.2.2 in [A\). We are grateful 
to A. Wakolbinger for proposing this interpretation. 

• Another interpretation is given in [19]. The Yaglom distribution of 
the process Xt, defined as the limit law as t oo of conditioned 
on > 0, is the exponential distribution with parameter 2\fi\/c (see 
Proposition 13.31 below, with ^ = 0). The Gamma distribution appears 
as the size-biased distribution of the Yaglom limit (P*(xoo G dr) = 
rP{Y e dr)/P{Y), where Y ~ £^xj»(— )), which is actually a general 
fact ([la Th.4.2(ii)(b)]). ' ^ 
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We have just proved that, for yU < 0, the law of Xt conditioned on Xt+g > 
converges to a Gamma distribution when taking first the hmit 6 oo and 
next the hmit t — oo. It is then natural to ask whether the order of the two 
limits can be exchanged: what happens if one first fix 6 and let t tend to 
infinity, and then let 6 increase? We obtain the following answer. 

Proposition 3.3 When fi < 0, conditionally on Xt+e > 0, Xt converges in 
distribution when t oo to the sum of two independent exponential r.v. with 
respective parameters — and ^(1 — e^^). 

Therefore, one can interchange both limits in time: 

lim lim P{xt G ■ | Xt+e > 0) = lim lim P{xt e ■ \ Xt+e > 0) = r(2, 

9^oo t^oo t^oo 9^oo C 



Proof First, observe that, by (1321) . 



A 2|/i| e'^* ^ X 
hm u. = and hm — - 



As in dH), it holds 

E(e-^"* I xt+e > 0) 



1 - Pixt+e = 0) 

, A + limT 

, A + liniT Un 

1 -exp(-x^Y3^jr{i+?)) 



Thus 

lim E(e-^"* I xt+e > 0) = ^r^e"'^^ lim el^l*K^+''"^^°° - wf) 

^ c ^-^e/ A + limx^^M^ A_ 

2H' Vi-|-(A + lim,_n^) 1 - ^^A 
_ 1 1 
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where the first (resp. the second) factor is equal to the Laplace transform of 
an exponential r.v. with parameter 2|/x|/c (resp. with parameter ^(1 — e^^)). 
This means that 

lim P{xt G ■ I Xt+e > 0) = Sxp{'^) ® Sxpi'^il - e^')). 

t— >oo C C 

It is now clear that 

lim lim P{xt e ■ I Xt+e > 0) = Sxp{'^) ® Sxpi'^) = 1(2, ^). 



-*oo t — ^oo 

Thus, the limits in time interchange. □ 



Remark 3.4 The previous computation is also possible in the critical case 
and gives a similar interchangeability result. More precisely, for any 6 > 0, 
Xt explodes conditionally on Xt+e > in P-probability when t +oo. In 
particular, for any M > 0, 



lim lim P{xt < M \ xt+e > 0) = lim lim V{xt < M \ xt+e > 0) = 0. 

(?— ►oo t— >oo t— >oo 6— >oo 



Remark 3.5 One can develop the same ideas as before when the branching 
mechanism with finite variance c is replaced by a /^-stable branching mecha- 
nism, < /? < 1, with infinite variance (see [2] Section 5 for a precise defini- 
tion). In this case, equation (IHTI) has to be replaced by = [iu\ — c{u\Y^^ 
which implies that 

^ _ A e'^* 

l + £^(l_e/^MO) 

Therefore, with a similar calculation as above, one can easily compute the 
Laplace transform of the limit conditional law of Xt when t ^ oo and prove 
the exchangeability of limits: 

lim lim E(e-^^* | Xt+e > 0) = lim lim PieT^''' \ xt+e > 0) ^ 



As before, this distribution can be interpreted as the size-biased Yaglom dis- 
tribution corresponding to the stable branching mechanism. This conditional 
limit law for the subcritical branching process has been obtained in |2Tj The- 
orem 4.2. We also refer to [19] Theorem 5.2 for a study of the critical stable 
branching process. <C> 
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3.2 Multitype irreducible case 

We now present the multitype generalization of Proposition l3. li on the asymp- 
totic behavior of the conditioned multitype Feller diffusion with irreducible 
mutation matrix D. 

Theorem 3.6 (a) In the critical case (fi = 0), when the mutation matrix D 
is irreducible, Xt explodes in P* -probability when t —>■ oo, i.e. 

Vi e {1, . . . , k}, MM > 0, lim P*{xt,i < M) = 0. 

t— >+oo ' 

(b) In the subcritical case < 0) when D is irreducible, the law of Xt 
converges in distribution under P* when t ^ oo to a non-trivial limit 
which does not depend on the initial condition x. 



Proof One obtains from 0241) that 



Dv^ - csupK^) v^<-^< Dv^ - cinfK^) ^;^ 

at t 



Then, by Lemma [2.41 

Therefore, in the critical case, vanishes for t large if A > 0, due to the 
divergence of miiU^^ds, which is itself a consequence of Lemma [231 (iii). 
If Aj = for some type i, by Lemma (23] (i) and the semigroup property of 
1 1—*- Ut, we can use once again Lemma [2.31 (iii) starting from a positive time, 
to prove that limi^oo = {i. Then, the explosion of Xt in P*-probability 
follows directly from f l23l) and from the fact that lim^^oo = 0. 

To prove (b), we study the convergence of v'^ := e~^*v^ when t oo. 
By ([3SD and Lemma 12.31 (ii) we know that t ^ is bounded and bounded 
away from 0. Fix e G (0, 1) and to such that J^^sup^M^j dt < e. Then, for 
any t > 0, 

< vi^, < e^'>-^^^% (36) 

and so, for any s, t > 0, 

- vlJ < sup I {e^'^I - e(^-^^)^) e^^->^'^% \ . 

<56{-l,l} 

By Perron-Frobenius' theorem, limj^oo e*-^"'^'^^* = P := {^if]j)ij and thus, 
when t — > oo, 

\vt,+t+s - Vto+t\ < (e^^-l)(5,^„,r/)e+|^i|o(l) 
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where the neghgible term o(l) does not depend on t^, s, e and A, since e <1 
and exp((D — is a bounded function of s. Therefore, {y^)t>o satisfies 

the Cauchy criterion and converges to a finite positive hmit v'^ when t oo. 

We just proved the convergence of the Laplace functional ( l23l) of Xt under 
P* when t — i> oo. In order to obtain the convergence in law of Xt-, we have 
to check the continuity of the limit for A = 0, but this is an immediate 
consequence of lim^^o lim^^oo = ^. 

Finally, letting t go to infinity in ( |36l) . we get 



\vt-Pvt\<i^-e-nPvto' 



where Pv^^ = {v^^,ti)C,. It follows that is proportional to ^, as limit of 
quantities proportional to C,. Therefore {x,v^)/{x,C,) = {v^,l) is indepen- 
dent of X and the limit law of Xt too. □ 

We can also generalize the exchange of limits of Proposition 13.31 to the 
multitype irreducible case. 

Theorem 3.7 In the subcritical case, conditionally on {xt+e, 1) > 0, Xt con- 
verges in distribution when t — > +oo to a non-trivial limit which depends only 
on 6. Furthermore, one can interchange both limits in t and 9 : 

lim lim P{xt e ■ \ {xt+e, 1) > 0) = lim lim P(xt G ■ | {xt+e, 1) > 0). 
Proof Following a similar computation as in the proof of Proposition 13.31 

lim E(e-(-'^) I 1) > 0) = lim exp(-(a:, .^)) - exp(-(x, .r^^- "^)) 

''t+e 



1 - exp(-(x, limx^^ g)) 



Since 



one gets 



Du^ - -{snpu') u^<-j^< Du^ - -(inf 



exp ( / supWg j (is )e^*A < -u^ < exp ( / infUg^ds 

^ 2 Jo i ' / \ 2 Jq i ' 



e^*A. 

This inequality, similar to flH^ . can be used exactly as in the proof of Theo- 
rem [32] (b) to prove that, for any e, there exists to large enough such that 
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and to deduce from flHTj) that := e~^*u^ converges as t ^ oo to a non-zero 
limit proportional to the vector ^. 

Moreover, because of ([5]), is increasing with respect to each coordinate 
of A. Therefore, it is elementary to check that t ^ limi^^ = sup„ m""*" 
is also solution of the non-linear differential system ([6]), but only defined 
on (0,00) (recall that, by Lemma fT^ (ii), lim^^^^w^ < 00 for any t > 0). 
Indeed, assume that bf = sup„ where d" = F{a") for a locally Lipschitz 
function F. Fix t such that bt < +00 and a small 1] > 0, and let F_ := 
inf |^_b^l<^ and F := sup|^_(,^|<^ There exists no such that, for n > 

no, — bt\ <Ti/2. Moreover, for any s in a neighborhood of t, |a" — 6j| < 77, 
where the neighborhood depends on F and F, but is uniform m n > uq. 
Therefore, for sufficiently small s and for n sufficiently large, F_ < {a^_^_g — 
o-t) / s < F . Letting n —>■ 00, s — > and finally rj ^ 0, since F — F — > when 
T] ^ 0, bt is differentiable at time t and bt = F{bt). 

Therefore, the semigroup property of the flow of (E]) implies that, for any 
t > 0, 

hm <, = nf-^-"^ 

SO that e~^^\im.)^^^u^ also converges as t ^ 00 to a positive limit tt^, 
proportional to ^ too. 



Hence, 

lim E(e~^^'''^^ | {xt+e, 1) > 0) = e 



= e-^^^^^^^^ (38) 

which is independent of the initial condition x. 

In order to prove the convergence in distribution as t — > 00 of a;t condi- 
tionally on {xt+e-i 1) > to a random variable with Laplace transform ( |38l) . it 
remains to prove the continuity of this expression as a function of A for A ^ 0. 
To this aim and also to prove the exchangeability of limits, we use the follow- 
ing Lemma, the proof of which is postponed at the end of the subsection. This 
lemma gives the main reason why the limits can be exchanged: is solution 
of the linearized equation of ([6]), and therefore, the gradient of with respect 
to A is solution of the same system of ODEs as v'^. The function was in- 
volved in the computation of limj lim^ F{xt G ■ | Xt+e > 0), whereas the gradi- 
ent of will be involved in the computation of limg lim^ P{xt G ■ | Xt+e > 0). 

Lemma 3.8 The function X ^ is differentiable, and its derivative in 
the direction rj, denoted by V^m^, is solution of the same differential system 
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as except for the initial condition given by V^Uq = rj. Furthermore, 
X 1-^ is dijferentiable too and its derivative in the direction rj, denoted by 
S/r^u'^ satisfies 

V^M^ = lim e'^^V^Ut. 

t— ►oo 

Since Vr,M^ satisfies the same differential equation as v^, in particular, 
V^M^ = and, with the notations of the proof of Theorem 13.61 V^m^ = v'^. 

It also follows from the above lemma that is continuous as a function 
of A. As a result, 



oo 



lim lim E(e-("*'^) | (x^+g, 1) > 0) = lim e'^^> ' 



Moo ',1) 

e ' 



smce 

lim e ^ M( * 



^ihn,^^«, = lim e-^^wf'"^-*-"* = lim e"'^* Jim uf^g = e^'w^. 



Finally, let us check that the limits in t and 6 can be exchanged. Since 
lim^^j^Mg ~ e^^u'^ = e^^{u'^, 1)^ when ^ oo, it follows from Lemma [XSl 
that 



lim e"^^ ^"- "-^^ = (V,«^, 1) = {vi, 1) = lim r(e-(-'^)), 

which completes the proof of Theorem 13.71 □ 



Proof of Lemma 13.81 The differentiability of with respect to A and the 
ODE satisfied by its derivatives are classical results on the regularity of the 
flow of ODEs (see e.g. Perko [26]). 

Moreover, since V^m^ and are both solution of the ODE ([21]) (with dif- 
ferent initial conditions), it is trivial to transport the properties of proved 
in the proof of Theorem 13.61 to VrjU^. In particular, e~^*V,,M^ converges as 
t +00 to a non-zero vector which is proportional to ^. We only have 
to check that Vj^tt^ exists and that = ^rjU'^- Moreover, as for ( l35l) . 

exp ( - I y supw^j ds^e^^rj < V^m^^ < exp -j^ J i^f ^s,, ds^e^^rj. 
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Therefore, since exp{Dt) > 0, 



which imphes that e~'^*V^M^ is uniformly bounded for t > 0, A > and rj in 
a compact subset of R''. 

Now, letting t +oo in (1H7|) one gets for any h > 0, 

^^^^ -ut-P r e~^'^V,uir''dr\ = \u'J'^ - - Pe-'^*°(<+^" - <)| 



< (1 - e-^^/2)p(g-/.io^A ^ g 



to ^ "to 



Letting £ — (and thus to ^ +oo) in the previous inequality, Lebesgue's 
convergence theorem yields 



h 



Therefore, is differentiable with respect to A and V^u^ = w^. The proof 
of Lemma 13.81 is completed. □ 



4 Some decomposable cases conditioned by 
different remote survivals 

In this section, we study some models for which the mutation matrix D is not 
irreducible: it is called 'reducible' or 'decomposable'. In this case, the general 
theory developed above does not apply. In contrast with the irreducible case, 
the asymptotic behavior of the MDW process and the MF diffusion depends 
on the type. 

Decomposable critical multitype pure branching processes (without motion 
and renormalization) were the subject of several works since the seventies. 
See e.g. [231 [IB [H EHl E31 EI] ■ 



4.1 A first critical model 

Our first example is a 2-types DW process with a reducible mutation matrix 
of the form 

^=("oo)' ^^^^ 

For this model type 1 (resp. type 2) is subcritical (resp. critical). Moreover 
mutations can occur from type 1 to type 2 but no mutations from type 2 to 
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type 1 are allowed. 

In this section we analyze not only the law P* of MDW process conditioned 
on the non-extinction of the whole population, but also the MDW process 
conditioned on the survival of each type separately. 

Theorem 4.1 Let P be the distribution of the MDW process X with muta- 
tion matrix flHUj) and non-zero initial condition m. Let us define P*, P* and 
P* for any t > and B e J^t by 



P*{B) = 


lim P(B 1 

e~*oo 




>o) 


P*{B) = 


lim P(B 1 




e,i,l)>0) 


P*{B) = 


lim P(B 1 


{Xt+e,2, 1) > 0) 



Then, all these limits exist and, for any t > 0, 



P*\ =p*\ =i^^p| (40) 



and PL =^^^^e"*PL. (41) 
(mi, 1) '-^^ 



Proof Let us first prove ( HTl) . Using the method leading to (fT2i) . we get 
that if mi ^ 



P*{B) = lim e(1b 



^Xf,limAi-^ooM^^''°^) 
(a;,limAi^oo u^+f^) 



The cumulant of the mass process satisfies for any A = (Ai, A2) 

= -"«M + "^t,2- 2K1) <l = ^l 

The second equation admits as solution 

= — (42) 
and can be computed explicitly if A2 = 0: 
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We then get 



P*(B) = e'''E(lB—), 



which yields (1411) . 

Concerning P*, remark first that it is well defined even if = (but 
mi 7^ 0) since particles of type 2 can be created by particles of type 1. 
We are going to prove (HOl) by a similar method as Theorem 12.21 Let us first 
compute ^. The matrix D has two eigenvalues, and —a, each of them with 
one-dimensional eigenspace. The (normalized) right and left eigenvectors of 
the greatest eigenvalue = are respectively ^ = (|; ^) and r] = (0; 2). Since 
^ > 0, the proof of Lemma [2.31 (and therefore Lemma [2.31 itself) is still valid 
for this specific matrix, except for the assertion (i), which has to be reduced 
to the following: if A2 > 0, then, for any t > 0, > (if A2 = 0, then 

Therefore, as in the proof of Theorem 12.21 we can prove that 
VA > 0, Vt > 0, - e^*A| < K\\X\\te^'X 
and thus that, for 6 and to such that \\u^J\ < 1/ K{9 + 1), 

(a,<+e) (a,e™0 



(&,e^(^+*)<) 

^ 2i^-||a||||<||g||e^Vj| 2i^-||a||||&||||e^Xl|||<||(t + g)||e^(^+^)< 



(6,e^(^+*)0 (6,e^(^+*)02 

(44) 



Let us compute explicitly the exponential of the matrix Dt. Since 

1 -1 




[-aYN where N = { \ \ \, then 



e^* = P + e-"*iV, with P = {^iVj)i<i,j<2 



1 
1 



This has the same form as fllSp in the irreducible case, except that P ^ 0. 
Because of this, we cannot obtain a bound for (14^ uniform in A G as 
in the proof of Theorem 12. 2[ However, we can restrict to a subset of for 
which the convergence is uniform and which covers the two limits involved 
in the computation of P* (Ai = and A2 +00) and P* (Ai +00 and 
A2 +00). 

This can be done as follows: if Ai = A2, then u^^ = u^2 for any t > 0. 
Therefore, for any A 7^ {0} such that A2 > Ai, u^2 > > for any t > 
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(if at some time these quantities are equal, they remain equal for any larger 
time). Then, since 

this quantity converges when 6* — > oo to (m^^ 2j'^to.2)y uniformly in A such that 
Aa > Ai > 0. 

From this follows as in the proof of Theorem 12.21 that 

lim 7 T — - = - — — 

uniformly for A in the set of (Ai; A2) 7^ such that A2 > Ai > 0. This ends 
the proof of (HO]). □ 

From the local density of P* (resp. P*) with respect to P, we easily obtain 
as in Section 12.21 the following expressions for the Laplace functionals of the 
different conditioned processes. 

Theorem 4.2 The probability measure P*(= P*) is characterized by 

V/ e a(P, P2)+, r(exp-(X„/) I Xo = m) = i^^l^^^e-^"^'^*^^ 

where Vtf is the unique semigroup solution of the PDF f l22p . The probability 
measure P* is characterized by 

V/ea(R,R')+, r(exp-(Xi,/) |Xo = m) = ^pi^e"V<'".^*^> 

(mi,l) 

where Vtf satisfies the same PDE as Vtf except for the initial condition 

K)/ = (i;0). 

4.2 Long time behaviors of the Feller diffusions 

Let us now analyze the long time behavior of the various conditioned MF 
diffusions. 

Proposition 4.3 (a) The first type vanishes in F* -probability whent — >■ 00, 
that is 

> 0, lim P*(xj,i > e) = 0. 

(b) Under P* , the first type converges in distribution to the probability mea- 

sure r(2, 2a/c) 

lim P*{xt,i e ■) = r(2,2a/c). 

(c) The second type Xt,2 explodes in P*- and in ¥* -probability when t —>■ 00. 
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with 



Proof We first compute the vector := V^A, A G R^. 

r (exp A) I xo = x) = i^e-(^'"^') 

= -tt^M + "^^2 - CMt^^M> ^0,1 = 1/2, 

^ = -CW',2<2, <2 = 1/2- 

Therefore, replacing by its value obtained in 

1 



2(1 + f A2t)2 

and 



In particular, if A2 = 0, 2^'°'' = 1/2 and one gets from the explicit 
sion (USD of 

(A„o) _ 1 cAi(l + |^AOte-"* + gAfe-^"^ 



2 2(l + ^Ai(l-e-0)' 

Similarly, for {)^ := VjA, A G R^, we get 



^^_yt)_^at-{x,ul) 



E*(exp -(Xi, A) I xo = x) = V^e"*e 



with 



-a^i^i + avl^ - cul^vl^, v^^^ = 1, 



4i 

= -CUt2Vt,2, <2 = 0. 



Therefore, -0^2 — 



t'^;^ = exp — at — c J u^ -^ ds 
In particular, if A2 = 0, using P3l) again. 



„— at 

^,(Ai;0) _ e 



(l + ^Ai(l-e--))' 
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Now (a) and (b) can be deduced from the facts that limt^ooVt'^^''^^ = ^ 
and hm^^oo e°*#^'°^ = (1/(1 + cXi/2af; 0). 

The explosion of Xt,2 in P*-probabihty in (c) is a consequence of the fact that 
hm^^oo = (exp(— c /q°° Ms^i c^s) ; O) = (0;0), by Lemma 

Finally, it follows from ( l45l) that 



m 



-at p-atjl 



e e ' a 



cA2(l + fA2t) 

where the last inequality is obtained by splitting the integral over the time 
interval [0,t] into the sum of the integrals over [0, |] and This implies 

the first part of (c). □ 

We interpret this proposition as follows. Conditionally on the survival 
of the whole population, the weakest type gets extinct and the strongest 
type has the same behavior as in the critical monotype case. Conversely, 
conditionally on the long time survival of the weakest type, the weakest type 
behaves at large time as in the monotype sub critical case and the strongest 
type explodes. 



4.3 A more general subcritical decomposable model 

We consider a generalization of the previous model. The mutation matrix is 
now given by 

«=("o_;) w 

where a > (as before) and /? > with [3 ^ a. 



In this case, the whole population is subcritical. Here again, mutations 
are only possible from type 1 to type 2. \i [3 < a, type 2 is "less subcritical" 
than type 1 (as in the previous case) but if a < /?, type 1 is "less subcritical" 
than type 2. We will see below that the behavior of the various conditioned 
processes is strongly related to the so-called dominating type, which is the 
first one if a < /3 and the second type ii (3 < a. 

Before treating separately both cases with different techniques, we define the 
common ingredients we need. 

We can easily compute the normalized right eigenvector ^ for the greatest 
eigenvalue /i. If /? < a, = — /3 and ^ = 2^^(«;« — /5) and if a < /?, 
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fi = —a and ^ = (1;0). We can also explicitly compute the exponential of 
the mutation matrix: 







The cumulant of the mass process satisfies 



dull 



-aut, + auto - ^Ki)^ <i = Ai 



du 



5 ' 2 ' (49) 



Thus M^_2 is given by 



One can compute explicitly only when A2 = 0, and in this case, as in 

= 7T aild Ul'o'^^ = 0. 

l + ^Ai(l-e-"*) 
We now consider the system of equations 

dht 1 , , ^ , 

which solutions are given by 

2/3 

and ^ 

i = e-'''-''^'»<^'^'{ho,i + a e"^+"^o <i'^"/i,,2c/s). (53) 

We denote as before by f^, or -y^ the respective solutions of (I5T1) with 
initial conditions = ^, {)o = (1; 0) and -Uq = (0; 1). 



(52) 
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4.3.1 Case p < a 

We now identify the laws obtained by conditioning with respect to the various 
remote survivals. 

Theorem 4.4 Let P* (resp. P*,P*j he the conditioned laws defined in The- 
orem WA\ where P is the law of the MDW process with mutation matrix given 
by fl48p with (3 < a and non-zero initial condition m. It holds 

P*| = p*| ^i^iiilg/^* p| 

\rt Ij^t (m, ^) 



and P*| =i^^e"*PL (^/mi ^ 0). 



Sketch of the proof The greatest eigenvalue of D is /i = — /5, the nor- 
malized right eigenvector for /i is ,^ = 2^;zrg(a; a — [3) and the normalized left 
eigenvector is = (0; ^^E^)- 

As in the proof of Theorem 14.11 ,^ > 0, so that Lemma 12.31 holds (except 
assertion (i) ) and we can use a similar method. The only difficulty is to 
find a domain C such that, for each initial condition X & E, the cu- 
mulant semigroup takes its values in E and {A1/A2, A e E} is bounded. 
To this aim, one can check that, if < u^i = -^^u^2 at some time t > 0, 

then ^ < Therefore, if < Ai < ^Aa with A2 > 0, one has 

< < -^^u^^2 any positive t. □ 

Let us now analyze the behavior for large t of the mass process under 
the three measures P*, P* and P*. Since {)o,2 = 0, -04^2 = as in Section W7I\ 
and ( H6l) holds. Therefore the behavior of Xt under P* is exactly the same as 
for /3 = 0, treated in Proposition 14.31 (b) and (c). 

The long time behavior of Xt under P* is different from Section 14.21 and is 
given in the following proposition. 

Proposition 4.5 (a) The first type vanishes in ¥* -probability when t 00. 

(b) Under P* , the second type converges in distribution to the probability 
measure r(2,2/3/c) 

lim P*(xi,2 G ■) ^=^r(2,2/?/c). 
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Proof Since 

we have to compute linii^oo f^e^*. 

For the proof of (a) we remark that, from (!52l) and ( !53l) . 2^'°'' = ^^2 e"''* and 



V- 



(Ai;0) _ ^-at 



Since 

ft 



one obtains 



,(Ai;0)^/3t _ {a-p)t \ 

2^' 



(l + ^A,(l-e--))^ / ^ 



The integral can be computed exphcitly and is equal, for t large, to 

a — p la 



Thus 



lim.a-V* = ^6 = ei. 
For the proof of (b), it suffices to show that 

From (l52l) . it is clear that lim^^oo "^^t 2^^^ e^* = n-i-J-\ 12 ^- It then remains to 
compute the limit, for A2 > 0, of v'^^^'^^ as t ^ 00. Using ( l53i) . we get 

■/o (l + ^(l-e-/^^))^ 
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The first term is 0{e '^^*) and goes to as t — > oo. Tlie limit of tfie 
integral can be computed as follows: 



{l + f X,(l-e-^^)y 



ds — 



2P' 







< K 



1 + ^ 

^ 2/3 



ds 



V 



1 - 



cA2„-/3t/2 _ 
2/3 ^ ^ 

1 + ^As 



where the positive constant K may vary from line to line and where the 
last inequality is obtained by splitting the integration over the time intervals 
[0,|]and[|,t]. 

Now, by Lemma [273] (ii). Imit^oo jij^ u^^^ds = 0. Therefore, 



lim v^^e'^^ 



"6 



t^oo 



-(a-/3)s^^ 



6 



(1 + ^^2) 



as required. 



□ 



Here again, this result can be interpreted as follows: for i = 1,2, condi- 
tionally on the survival of the type i, this type i behaves as if it was alone, and 
the other type j explodes or goes extinct according to whether it is stronger 
or weaker. 



4.3.2 Case a < (3 

When a < f3, the greatest eigenvalue of -D is /i = —a and the normalized 
right eigenvector to /i is ^ = (1; 0). In particular, ^ / 0, so that Lemma [2.31 
does not hold and we cannot use the previous method anymore. However, 
in our specific example, can be explicitly computed, and, by (H9l) . u^^ is 
solution of the one-dimensional differential equation 

dyt c 2 , aA2e-^* 

^ = - 2^* + l + ^A2(l-e-/^r ^ ^ 

This equation can be (formally) extended to the case A2 = 00 as 

dyt c 2 , 2a/5e-^* 

IF = ' 2^* + cd-e-z^n - ^''^ 
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The following technical lemma gives (non-explicit) long-time estimates of the 
solutions of fl3^ that are sufficient to compute the various conditioned laws 
of the MDW process. We postpone its proof at the end of the subsection. 

Lemma 4.6 For any A2 G [0, 00], let 3^(A2) denote the set of solutions yt 
of fl3^ (or of fl3^ if ^2 = co) defined (at least) on (0, 00). For any y G 
3^(A2), the limit C{X2,y) '■= lim^^oo l/tc"* exists and satisfies 

0< inf C{X2,y)< sup C(A2,y)<+oo 
A2>i,yey(A2) A2>i,j/ey(A2) 

We now identify the laws obtained by conditioning P with respect to the 
various remote survivals. 

Theorem 4.7 Let P* (resp. P*,P*) be the conditioned laws defined in The- 
orem where P is the MDW process with mutation matrix given by ( HHl) 
with a < (3 and initial condition m = {mi, 7712) with nii 7^ 0. It holds 

p* I _ p* I _ p* I _ '^^ gOt p I 

\Tt \Tt \Tt (m,^) 
When mi = and m2 7^ 0, 

p*| ^ p*| _ {^t,2, 1) ^(3t p| 

and P* is not defined. 

Proof Our usual method consists in computing the following limits as Q 

goes to +00 : 

illllAi^oo^e ilHlAa^ooMg Wm.\^^\^_,aoUQ ^^^^ 



(x,lim.,_«;if)' (x,lim.,_«S^V (x,lim.,.,_«;^f^) 



It is elementary to prove that, as monotone limits of solutions of (!54|) . the 

l^Ai— >oo 



function Q 1— >• limA,_^oo ^^i^i^'^^ is still solution of fl54l) with A2 = 0, and the 



functions Q ^— limAj^oo j ^ and 6* 1— > limAi,A2^oo 1 ' ^ are solutions of (l55l) 
(a priori only defined for t > 0). 

Therefore, we can use Lemma 14.61 and the explicit formula flSUI) for 
compute the three limits of fl56l) . In each case, the dominant term is the one 
including u\y, except when m\ = 0, where the only remaining term is the 
one including u^2- ^ 

Finally, we give the long time behavior of the mass process under P* 
(which is equal to P* and P* when this last measure exists). 
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Proposition 4.8 (a) If mi ^ 0, the laws under P* of the mass of both types 
Xt^i and Xt^2 converge when t ^ oo. More precisely 

lim ?*{xt,i e ■) = r(2,2«/c). 

Xt,2 converges in ¥* -distribution to a non-trivial (and non-explicit) dis- 
tribution on [R+. 

(b) //mi = (m2 ^ 0), Xt,i = P* - a.s. and 

lim P*(xi,2e-) = r(2,2/3/c). 

t— »oo 

Proof With the previous notation, when mi ^ 0, 



E*(exp-(a;i,A)) = )3^e"V 



{x,vl^ 

Since Vt^2 = and Vt^i = exp(—at — c J^Ug ids), 



oo 



lim e"*'Of = ( exp — / u^ids; 

t^oo V Jo ' 



where exp — u^ids G (0, 1) by Lemma [4. 61 In order to prove the conver- 
gence in distribution of Xt,2, it remains to prove that 



oo 



lim / u?[^''^ds = 0. 
^2^0 Jo ' 

Because of (EH), Uti — "(^'^ti + «A2e~^*. Therefore, it is easy to check that 
«t^i <iP + 2)A2e-^* for all t > if Ai < (/5 + 2)A2 (simply differentiate the 
difference). This implies the required result. 

When A2 = 0, the computations can be made explicitly as in the proof 
of Proposition 14.31 (b) and give the usual Gamma limit distribution for xt^i 
under P* when t +00. 

If mi = (i.e. Xi = {mi, 1) = 0), 

r(exp-(xi. A)) = i^e^*e-(^'^') = t;,\e^*e-(^'"* ^ 
[x, i)^) 

and the computations are the same as in the monotype case. □ 

We interpret this last result as follows: when the first type is present, it 
dominates the asymptotic behavior of both types, since its subcriticality is 
weaker than the one of the second type, although mutations from type 2 to 
type 1 do not occur. 
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Remark 4.9 Theorem 14.71 and Proposition 14.81 are still valid in the re- 
maining case a = (3. In this case, it is actually possible to be more pre- 
cise than in Lemma 14.61 by proving, using a similar method, that a so- 
lution yt to ( 15^ or ( l55l) (with a = [3) satisfies yt ~ C{\2)te~"^^ where 
C(A2) = a\2/{l + cA2/2a) if A2 < oo and C(oo) = 20^/0. 

Proof of lemma 14.61 Let Zt := e°*?/j. It solves the equation 

dzt c 2 , aAse-^'^-")* 



—.e-'4+ . , n (57) 

2/3' 



dt 2 * I + S^X^il - e-P' 



Let us first prove that Zt is bounded for t G [1, cxd), uniformly in A2 and 
independently of the choice of the solution yt of or (1331) defined on (0, cxd). 
For any t > 1 and A2 G [0, cx)], 

dz^ ^ 2a(3 ^_(i3_a)t 
dt - c(l - e-^) 

Since the integral of the above r.h.s. over [1, 00) is finite, we only have to prove 
that zi is bounded uniformly in A2 > and independently of the choice of 
yt- Now, for any t G [|, 1] and A2 G [0, 00], 



dzt 2aP e-^/^-")* ce"" , 
< z ^ —z: 



dt - c l-e-/3/2 



In particular, for any t G [|, 1], the first term in the r.h.s. above is bounded 
and bounded away from 0. Thus, there exists a constant K such that, if 
Zt > K and | < t < 1, ^ < —^^^z^. Therefore, distinguishing between 
zi/2 ^ K and 2:1/2 > one obtains 

-21/2 4e" 

zi < yK< y K =: K' < 00. 

- 1 + fe— 2:1/2 - c 

Second, it follows from ( 1571) and from the boundedness of zt that 
\dzt/dt\ < if"(e~^^~°^* + e~"*) for some constant K" for any t > 1. Therefore, 
Zt converges as t — 00. Moreover, this function is uniformly bounded from 
above for t > 1 by some constant K'" independent of the particular function 
Zt considered. Therefore, since dzt/dt > —cK"'e~°'^Zt/2 for t > 1, the limit 
of Zt when t — +00 is also greater than 2:1 exp(— cK"'e~'^/2). 

Then, it only remains to prove that Zi is bounded away from 0, uniformly 
in A2 G [1, 00]. This follows from the fact that, for any A2 > 1, there exists a 
constant M > independent of A2 and t such that, for t G [|, 1], 

dzt , , c 
— ^ > M - -zf. 
dt - 2 * 
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This imphes that there exists M' such that, if Zt < M' for t G [|, 1], dzt/dt > 
M/2, and thus 

zi > {zi/2 + — ) A M' > — A M' > 0, 
which completes the proof of Lemma 14.61 □ 



4.4 Exchange of long time limits 

As in the irreducible case, one can study the interchangeability of the long 
time limit {t +oo) of the conditioned Feller diffusion and the limit of long 
time survival {6 —>■ +oo). The same method as in Proposition 13.71 yields, for 
i,j G {1,2} and A G IR+, 

hm E(e-^-''» I xt+e, > 0) = lim ^-^^^ 

t^oo t^oo (a;, lim^_,^ U^^g) 

where = (A; 0) and A^ = (0; A). 

However, the computation of these quantities requires precise informa- 
tion about the behavior of function of its initial condition A for t 
large. The cases we could handle are the one with non degenerate limits. 
In the model introduced in section 14.11 it corresponds to i = j = 1 and the 
computation reduces to the monotype case studied in Proposition 13.31 In 
the model introduced in section 14.31 with P < a, it corresponds to z = j = 1 
and i = j = 2, and with a < (3, to all cases. For i = j = 1 the proof is 
based on explicit expressions like in the monotype case, and for the other 
cases, the arguments are similar to those used in the proof of Theorem 13.71 
(except in the case a < P and mi = 0, where the computation can also be 
done explicitly). To summarize: 

Proposition 4.10 In the cases described above, one can interchange both 
limits in time: 

lim lim P{xt^i G ■ | Xt+e^i > 0) = lim lim P(a;(,i G ■ | Xt+e.i > 0) 

0— »oo t— >oo ' ' t— »oo (?— >oo ' 

= r(2,2a/c) if0<l3<a, 
hm hm P{xt,2 G ■ | Xt+e,2 > 0) = lim lim P{xt,2 G • | Xt+e,2 > 0) 

r(2,2/?/c) ifO<f3<a 

or < a < P and mi = 0, 
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lim lim P{xt,i G • | Xt+e,i > 0) = lim lim P{xt,i G • | Xt+e,2 > 0) 

p— >oo t— »oo 6— >oo t—^oo 

= lim lim P{xt^i G • | Xt+0,1 > 0) = lim lim P{xt,i G • | Xt+0,2 > 0) 

t— >oo 9— >oo ' ' t^oo 9^00 ' ' 



= r(2,2a/c) z/0<a</3 



anc? 



lim lim P{xt,2 G • | Xt+e,i > 0) = lim lim P{xt,2 G ■ | Xt+0,2 > 0) 

d^oo t— »oo 6— >oo t— >oo 

= lim lim P(xt,2 G • | Xt+e,i > 0) = lim lim P{xt,2 G • | Xt+0,2 > 0) 

t— »oo (?— >oo t— >oo 0— >oo 

if < a < (3 and mi ^ (in this last case the limit is not known explicitly). 
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